Has the nonlinear Meissner effect been observed? 
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We examine recent high-precision experimental data on the magnetic field, H, dependence of the 
penetration depth \(H ) in YBa2Cu307_i (YBCO) for several field directions in the a — b plane. In 
a new theoretical analysis that incorporates the effects of orthorhombic symmetry, we show that the 
data at sufficiently high magnetic fields and low temperatures are in quantitative agreement with 
the theoretical predictions of the nonlinear Meissner effect. 



72.40.Hi,74.25.Nf,74.20.De 



It is widely accepted jjj that the symmetry of the order parameter (OP) in high temperature superconductors 
(HTSC's) is at least predominantly d-wavc, vanishing at nodal lines approximately ninety degrees apart in a quasi 
two-dimensional Fermi surface (FS). Many details of the OP in these materials remain quite unclear, however. Are the 
^ nodes exactly at right angles? Are they true nodes or only very deep minima, "quasinodes" ||]? Addressing these and 
similar questions is important for obtaining clues about the nature of the superconductivity. They are particularly 
difficult to answer for the bulk OP, which may well differ || from the more easily observed surface state. To probe the 
bulk OP it is best to use electromagnetic techniques, since electromagnetic fields penetrate the sample over a depth 
of the order of the London penetration length A which is, in these materials, several orders of magnitude larger than 
the coherence length £ which characterizes the range of typical surface probes. Indeed, one of the early key results 
Q in support of bulk ei-wave superconductivity was the measurement of the linear temperature dependence of the 
^ ' penetration depth in a high-purity YBCO single crystal. However, such data indicated only the existence of nodal 
lines without the angular resolution needed to identify their position. Consequently, intensive efforts to precisely 
O ' determine the structure of the bulk OP have continued. 

,<->: it was first pointed out || years ago that nodes in the OP yield distinctive and measurable nonlinear effects in 
the field and angular dependence of the penetration depth when the superconductor is in the Meissner state. In 
t— I , subsequent theoretical work ^-|| more emphasis was placed on the existence, due to this nonlinear Meissner effect 
(NLME), of a component of the diamagnetic moment normal to the applied field, and on the torque associated with 
this transverse component. These phenomena were deemed to be easier to measure than the changes in A itself. It 
was shown || in this context that the NLME can be used to perform node spectroscopy, that is, not just to infer 
the existence of nodes, but to locate their positions on the FS, and to determine whether they are true nodes or not. 
Thus, the NLME is potentially a very important tool for the study of the pairing state in HTSC 's, as well as other 
materials in the ever increasing list of those for which the proposed OP leads to an energy gap with nodes. Yet, 
the experimental situation is rather confusing. The best experimental effort to measure the transverse diamagnetic 
moment [[To| in YBCO was inconclusive. Subsequently, results |ll[ for the magnetic field dependent penetration depth 
to a precision |l2| of ~ O.lA were reported. The NLME should be observable in such a high precision experiment, 
more precise than existing transverse moment measurements. Unfortunately, no theory of the NLME contribution to 
the penetration depth for orthorhomobic structures such as YBCO was available when Ref. |tl| was written. Only 
very recently [p"3f have the necessary calculations been performed. This has resulted in contradictory claims as to 
whether observed results are in agreement or not with NLME theory. Thus, a certain amount of skepticism has 
developed as to the observability of the NLME. 

In this paper we show that measurements of the field dependent penetration depth A(V>, H ) as a function of the angle 
ip that an applied field H in the a — b plane forms with the a-axis, must be analyzed very carefully. The anisotropy of 
the linear penetration depth tensor has a drastic effect jl3| on the NLME for X(ip, H). This may have been overlooked 
because the anisotropy effects in the transverse moment are known j^] to be relatively minor. One must also eliminate 
several other factors that may mask the signal at low fields and which are very difficult to account for theoretically. 
Thus we reanalyze here the best data available for the penetration depth in YBCO. We find that, although some 
questions remain, the low temperature data are quantitatively in agreement with theoretical expectations for the 
NLME in this material. 

We focus here on YBCO, the most experimentally studied [ p~0|j3~T| J]~^ , [3~5[ ] HTSC in this context. Hence, the relevant 
material parameters are well known, thus reducing the uncertainty in the fitting procedures. We perform our analysis 
primarily on the most complete available high resolution data of Ref. [pd| which includes results for four different 
directions of the applied field in the a — b plane. 
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The angular and field dependent increase in the penetration depth due to the NLME for materials with orthorhombic 
anisotropy of the YBCO type was first calculated in Ref. [O. The details will not be repeated here. The sample is 
assumed to have its larger faces parallel to the a — b plane (this is the case for crystals grown by the usual methods) 
and thickness large compared with the penetration depth. One has for the quantity AX(ip, H) = \(ip, H) — X(ip, 0): 

AA(V,H) = i-^W). (1) 
d n 

Here A is the geometric mean of the two in-plane principal values, A Q and A;,, of the zero field penetration depth 
tensor, Hq is a characteristic field of order <I>o/7r 2 A£ (<&o is the flux quantum), and y carries the angular dependence. 
The orthorhombicity, very important in this case, is incorporated into y through two parameters [ fl3| |: one is the 
ratio A = A a /A&, and the other is the angle a that the Fermi velocity at the node located in the first quadrant forms 
with the a direction. Because of the orthorhombic distortion of the FS, this angle does not have to exactly equal 7r/4 
even for a pure d x i_ y i state, while the quantity A, for YBCO, considerably exceeds unity. Here we take the zero field 
quantities A a and A& fixed at their experimental jl6| values (1050A and 1575A), giving A = 1.5. In this case, the full 
expressions jl3],[l7| for y(ip) simplify somewhat and can be written as: 
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where the angle tpi is given by ipi = arctan(tana/A). 

Because of the orthorhombicity, the angular dependence of A\(ip,H) is quite different [|l3| for A = 1.5 than that 
found for the tetragonal case (A = 1, a = 7r/4). This is unlike the situation for the transverse magnetic moment [fj|, 
which from symmetry considerations vanishes when H is along the a or the b axis. Thus, moderate orthorhombic 
anisotropy induces only a relatively minor distortion in the curve between ip = and ip = n/2 since these points are 
so to speak, anchored. This is not the case for A\(ip,H): when the field is applied along -0 = the currents flow 
over a region of thickness determined by A& while if ip = 7r/2 the relevant skin depth is A a . The effect is nonzero, and 
different, in either case. This difference is compounded by the nonlinearity and the apparent overall symmetry of y is 
7r rather than tt/2 even for moderate orthorhombicity. Failure to take this into account leads to erroneous conclusions 
concerning the angular dependence of A\(ip, H). 

To analyze data in terms of Eqs. (^) and (||), additional considerations are needed. These expressions, indicating 
that AA is proportional to H, are valid at low temperature. Here, "low" temperature is a field dependent concept. The 
characteristic temperature separating the high and low T regimes is [||J7) T*(H) w A (H/ H ), where A is the gap 
amplitude. At any finite T, the validity of the above equations will break down at sufficiently small H. Further, the 
effect of impurities is not included. For the clean samples used in experiments |l(],[ll],[l4j this should affect [^|J^] only 
the small field results. The same is true of possible nonlocal effects |lj|. If they are present at all in this geometry [fl9|| , 
they would affect results at fields below [^(J 20 gauss. Ideally, one would like to take into account all of these effects 
by modifying the above formulas. However, it is not feasible at present to include all of these factors simultaneously 
in a reliable manner. It is therefore best to perform the analysis in a consistent manner in terms of data in the higher 
range of fields available, where these additional effects are all weak. 

In Fig. ^] we show best straight line fits to the 1.2 K data of Ref. |llj for H along the a and b directions. All 
data in the range H > 60 gauss are included in the fit. The cutoff of 60 gauss was chosen as the point below which 
deviations from a straight line begin and it will be shown below to lead to a consistent interpretation. The straight 
line does not intercept the origin of the original plot, which has to be shifted downward. This is as expected, since the 
experimental AA includes the previously mentioned temperature pi] ], impurity, and possible nonlocal effects which 
increase this quantity with respect to the theoretical, clean, zero temperature, local value. The shift is small, of order 
1 A, indicating that the sample is clean and any such spurious effects are small. The two slopes of the lines obtained 
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from these fits are the quantities y(0)X/H o and y(ir /2)A/ 'H respectively, where from Eq.(||), y(0) — (2/A 2 )sin 3 a, 
and y(n/2) = 2A 2 cos 3 a. Thus, the ratio of the two slopes depends only on a since we have used the independent 
experimental value A = 1.5. We then determine the a that fits this ratio and subsequently find the characteristic 
field Ho from either one of the slopes. The results are very sensible: we obtain a = ir/4 + tt/17 and Hq = 5660 gauss. 
The value for the angle between the Fermi velocity at the node and the a-axis exceeds 7r/4 by a small amount, as 
one would expect for a pure d x 2_ y 2 pairing state and a tight binding FS with a slight orthorhombic distortion. The 
value of the characteristic field is consistent with expectations ^|J9) and also with our cutoff choice for the field: in 
the range of fitting we have H/Hq > 0.01. This means that in this field range, the characteristic temperature T*(H) 
introduced above is of order of 4 to 12 K. Hence the 1.2 K data included in the fitting are in the low temperature 
regime and the procedure is consistent. 

Up to now, we have, however, fit two quantities with two parameters, although the reasonable values obtained for 
these parameters are encouraging. To go beyond, we now use the obtained values of a and H to plot the predicted 
slope of the high field data at \tp\ = tt/4 without any additional parameters. This is done in Fig. |[ Experimental 
results for fields applied in the ip = ±7r/4 directions are included. These results ought to be identical (even with 
the orthorhombic distortion) and their small discrepancy reflects systematic errors in the experiment. Nevertheless, 
the fit is excellent in the high field range. We also plot, (inset) with these parameter values, the predicted angular 
dependence of AA for YBCO. One can see that, with the orthorhombicity, this angular dependence differs considerably 
from that obtained for a tetragonal system, also plotted for comparison. The actual curve is not symmetric about 
7r/4 and its maximum is much less pronounced than that for the tetragonal case, which is characterized Q by a factor 
of \/2 between maxima and minima. Because of these differences, the attempt made by the authors of Ref. jll| to 
reconcile the angular dependence of their data with the theory of the NLME in a tetragonal system, had to fail. 

In Fig. ||we compare the theoretical results with other more recent data |l4| on YBCO at T=1.4 K. The parameters 
used are exactly the same as previously obtained. No new fits were performed. Results for the two directions available 
(field applied along the two principal axes) are shown. This data is in a more restricted, lower field range, and it has 
considerably more scatter than that of Bidinosti et al |ll[ . All that can be said with certainty is that it is consistent 
with the NLME theory with the same parameter values. 

In summary, the main result of the analysis presented here is that the best low temperature, high field, data p"l| ] 
on the nonlinear penetration depth in YBCO is in quantitative agreement, in its magnitude and angular and and 
field dependence, with the NLME theoretical expectations. Other data pj are also consistent with theory. Failure to 
observe the NLME in the transverse moment |l0| seems to be attributable to the actual precision in that experiment 
being just slightly less than what was in fact required. 

Two remarks must be added: first, the crucial influence of the orthorhombic anisotropy in the angular dependence 
of AA, which becomes very different (see Fig. |^) from that found for tetragonal symmetry, must be emphasized. 
Second, one sees the need to finesse the temperature, impurity and possibly other problems associated with smaller 
fields by obtaining and analyzing data at the hig hest possible fields below that of first flux penetration. Fortunately, 
this field is in the range 200-400 gauss 0,|ll],|§ for typical YBCO crystals. 

The question of the temperature dependence of the results [fi"l| , [l4t is less clear and needs further discussion: results 
obtained at 7 K for the same sample mainly discussed here are flll|| not substantially different from those at 1.7 K. 
With the characteristic temperature T*(H) in the range estimated above, it can be that the high field results are not 
yet affected by the temperature at 7 K while those at low fields are dominated by largely temperature-independent 
impurity effects. Indeed, it appears that a straight line fit to the 7K data at the highest fields (see Fig. 4 of Ref. [ Q| ) 
has a larger (in absolute value) vertical axis intercept than that for the 1.7 K data, which would be consistent with 
this scenario. Nevertheless, the weak temperature dependence of the data will remain a puzzle so long as a rigorous 
calculation including impurities, temperature, and possibly, nonlocal and other effects is not feasible. It is possible 
that these effects combine to yield a temperature dependence weaker than what the naive theory would predict. 

Finally, our analysis indicates that there is no significant is admixture to the d x 2_ y 2 gap, since such an admixture 
would lead to quasinodes and to ]13|| a considerable reduction in AA. Furthermore, the nearness of a to 7r/4 is 
consistent with the absence of a real s component as well. 

It would be desirable to perform measurements of AA in YBCO at additional values of the angle ip to verify in 
more detail if the curve in the inset of Fig. |^ is indeed closely followed. The comments made here on the proper 
way to analyze experimental data should also be taken into account in any attempts to use the NLME to elucidate 
the pairing states of other suspected unconventional superconducting materials for which it is estimated JDs] that the 
sensitivity of present penetration depth measurements is sufficient to probe the NLME. 

We thank CP. Bidinosti and A. Bhattacharya for very useful information about their respective measurements. 
This work was supported in part by Petroleum Research Fund, administered by the ACS (at Minnesota) and by 
DARPA and ONR (N00140010028) at Maryland. 
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FIG. 1. Magnetic field (H) dependence of AA (see text). The straight lines are fits to the 1.2 K data (circles and squares) 
of Fig. 3, Ref. |TTJ , for H > 60 gauss. Top: H applied along the 6— axis. Bottom: H along the a— axis. 

FIG. 2. AA(-ff) for H along ij) = |7r/4|. The straight line is the theoretical result at higher fields with the parameters 
extracted from Fig. |l|. Diamonds and triangles are the experimental data with H applied at tp = ±7r/4. Inset: predicted 
angular dependence (thin curve) of AA(i/j) including anisotropy. Bold curve: result for a tetragonal system. The amplitudes of 
both curves correspond to H = 180 gauss. 

FIG. 3. AA as a function of H. The straight lines are the theoretical results with the same parameters found in Fig. |IJ The 
symbols are experimental data of Ref. H|. Top panel: H at tjj = ±7r/2. Bottom panel: H along tp = 0, 7T. 
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Fig.l Halterman et al 




Fig. 2 Halterman et al 




Fig. 3 Halterman et al 




